Two photon interference and optical free induction decay 
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The two photon interference phenomenon is theoretically investigated for the general situations 
with an arbitrary input two photon state with and without photon polarization. For the case 
without polarization, the necessary-sufficient condition for the destructive interference of coincidence 
counting is given as the symmetric pairing of photons in the light pulses. For both case it is 
shown that the "dip" in coincidence curve can be understood in terms of the free induction decay 
mechanism. This observation predicts the destructive interference phenomenon to occur even for 
certain cases with separable input two photon state, but it can only be explained in terms of "the 
two photon (not two photons )interference ". 



PACS number: 05.30.-d,03.65-w,32.80-t,42.50-p 

In last two decades one of the most important pro- 
gresses in quantum optics is the experimental demon- 
stration of the Einstern-Podolsky-Rosen (EPR) [1] effect 
with the entangled photons. The creation of EPR en- 
tangled photons in a spontaneous parametric down con- 
version (SPDC [2]) not only implements the best test 
of Bell's inequality [3] with fewer loopholes [4-6], but 
also lays foundations for some quantum information tech- 
niques [7]. This important progress has inspired people 
to reconsider the profound observation that "...photon... 
only interferes with itself. " stated by Dirac in his famous 
book "T/ie Principle of Quantum Mechanics "[8]. 

Some interesting experiments of typical two pho- 
ton interferometer seem to illustrate the existence 
of interference between two different photons since a 
curve of "dip" was observed in the rate of two photon 
coincidence [9- 11]. However, a series of refined experimen- 
tal setups [12-14] concluded that Dirac is correct. They 
argue that "a two photon (not two photons) can also 
only interfere with itself" [14] to account for the exotic 
interference phenomenon. In this argument the crucial 
conception is the two photon or bi-photon, the insepa- 
rable photon pair depicted by a EPR state. In fact, the 
naive idea of "destructive interference between the idler 
and signal photons" can not give a correct prediction for 
a two photon interference experiment "with three arms" . 

If one believe that the idea of bi-photon is indeed neces- 
sary for all the two photon interference experiments, then 
two natural questions follow immediately : 1. Does there 
exist the separable pair of photons to give the destructive 
interference in the two photon coincidence counting rate? 
In this case the separable pair of photons can be imple- 
mented experimentally in the two independent pulses of 
photons. 2. If there is indeed a destructive interference 
for two independent pulses of photons, does it imply "in- 
terference between the idler and signal photons" ? 



I. TWO PHOTON INTERFERENCE OF TYPE I 

In this paper the two questions will be answered in an 
universal framework by considering a general input two 
photon state 



I*) = 



duj 



dCjf{io,Co)a+{uj)a+{^)\0). (1) 



The schematic setup for this consideration is illustrated 
in figure 1: 
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FIG. 1. Schematic of the two-photon interferometer 

the two lights coming from points a and b are mixed at 
point O on the 50-50 beam splitter BS and detected by 
the two photon counters at Di and D2 ■ Correspondingly, 
aj^(a;) ( a'j^{uj) ) is the creation operator for the photon 
of frequency lo { lu ) in the optical path mode A {B) from 
the point a (6) to O; fiuj^Co) is the distribution function. 
The optical path modes can be regarded as the modes of 
idler and signal lights in the usual SPDC experiment with 
a special distribution function f{uj,oj) oc d{uj + uii — luq) 
for a finite real number luq. The separable pair of photons 
corresponds to the case with the factorized distribution 
function f{uJ,u}) = fi{uj)f2{uj). 
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FIG. 2. The destructive interference 

With this configuration of gedankenexperiment, the 
above two questions can be asked in an unique way : 
what kind of two photon state (or what kind of distribu- 
tion function /(w,a))) can result in the curve of "dip" of 
the rate of two photon coincidence, which is illustrated 
in figures 2 as the typical destructive interference phe- 
nomenon. Here, D — li — si is the optical path differ- 
ence between aO and bO; aO = Si,bO — /i, DiO = S2, 
D2O = I2 are the lengths of optical paths. The rate Rc 
of coincident detection is measured as a function of the 
position oi BS and there is a "null" in coincidence i.e., 
Rc — a.t D — 0. This indicates a destructive interfer- 
ence. To explain the anti-exponential decaying behavior 
beside the "dip" (or anti-peak) point, we will resort to 
the free induction decay mechanism [15]. 

In the following discussion we do not make assump- 
tions about the light source, which may be the entan- 
gled photon pair created by the BBO nonlinear crystal 
or any two independent photon input pulse. We denote 
the positive frequency parts of the electric field at de- 
tectors Di and D2 by and -Ej"- We will choose a 
proper position as the coordinate origin so that we could 
compute the coincidence rate conveniently. According to 
Glauber's coherence theory [16,17], the probability per 
unit (time)^ that one photon is recorded at Di at time 
ti and another at D2 at time t2 is just the second order 
coherence function. 

a^\h,t2) = {^\E^(t2)E^{h)E+{h)E+(t2m (2) 

We recall that, when one consider the first order coher- 
ence the superposition of two obvious "paths" is avail- 
able. To describe the two particle interference phe- 
nomenon due to the second order coherence in a sim- 
ilar way, the generalized "path" is introduced by con- 
sidering Gt^l(ti,t2) = \'^{ti,t2)\^ where the bi-photon 
wave packet ^{ti,t2) = {0\E^{ti)E}{t2)\'^) is invoked 
as a two photon effective wave- function [14]. Most re- 
cently this result was generalized for higher order co- 
herence in time-domain [18]. Then, the two photon co- 
incidence rate can be given by the two time integral 

Rc = C^dhdt2\-^ih,t2)\^. 



In present discussion, in terms of the positive frequency 
part of the A (B) mode electric field at point a (6) with 
the spectral density g{(jj) 

/"OO 

E+,){T)=j^ du;g{u;)aAiB)He-'-' (3) 

the local field operators 

Et{h) = ^[^E:{n)+E+{r[)] (4) 

E+it2) = ^[Et{r!,)+^E+{r2)] 

are implemented by the 50 — 50 beam splitter with 
Sl+ S2 ^ I1+ h 

Tl = tl ,T2 —12 

C C 

, , h+ 32 , , S1+I2 

Tl = II = l2 

C C 

where c is the velocity of light. Then the two-photon 
wave function can be written explicitly as 

/■OO 

* = / F ■ [e~'^'^^^+'^^i'> ~ e-''^'^^'+'^^^^]dLjdCj (5) 
Jo 

with 

F = F{uj, = ^giu;)f{u, Cb)giCo). (6) 



II. THE NECESSARY-SUFFICIENT CONDITION 
FOR THE DESTRUCTIVE COINCIDENCE 

Now let's prove that the necessary-sufficient condition 
for the destructive coincidence i?c = at Z? = is 
f{uj,uj) = f{uj,uj). To show the sufficience we decom- 
pose the two-photon wave packet = -I- into two 
parts 

^1= dtjdwF-e-*("^^+'^^^)[l-e-*("-'^)*] (7) 

/•OO 

^-2 = / diudojG{oj,oj) ■ e-'('^^2+'^^0 
Jo 

where G{ll!, uj) — F{lu, lu) — F{lu, to) in the derivation, we 
have considered ti — t{ = — r2 = It is easy to 
see that *2 = and = v^iwhen G{lu,lu) = 0. On the 
other hand, 1'i1d=o = 0. Then, at D = 0, i?c = 0. To 
show the necessity, wc consider ]^']^ = l^ij^ + |'I'2|^+ 
2i?e*t«'2. Since ^-i = at Z) = , \-^\^\d=o = |«'2|^ 
So if Rc\d=o = 0, then \^\^\d=o = |*2|2 = i.e., *2 = 
0. Because G{uj,uj) is the Fourier component of 'i'2, the 
vanishing '^2 must lead to G = or f{uj,uj) — 

In fact this necessary-sufficient condition proved above 
has been implied in several discussions [19-22]. 
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To illustrate the main idea of physics implied by 
the above proof, we consider the simplest symmetric 
state I*) = \Sin,il)) = + ,which is 

an entangled state in frequency domain. Here, |0,f2) = 
a\{Q:)a^{Q)\0). The corresponding two photon wave 

f^^^) g (Q) g (^^ij has two 



packet = ^(0n 

parts </.^^ = er<''<+''^0 



and (j,^^ times 
the state density of the light field. Notice that the 
part (pQQ is contributed by the component |f2, f2) and 
(pQQ hy A straightforward calculation shows that 

^nn = at D = 0. We think this phenomenon can not 
be understood as "interference between two photons". 
The reason is, at _D = 0, neither 4>q{i nor (f)^-^^-^ is zero, 
but the sum of them is zero. In fact, the single term (pQQ 
( ) corresponds to the scattering of two independent 
photons in il) ( on the beam splitter. Due 

to their different frequencies fl and fl, two independent 
photons can not interfere with each other. In this sense 
the destructive interference can only be attributed to the 
fact that " the two-photon interferes with itself , which 
is similar to Dirac's statement for single photon. 

Since any symmetric two photon state j^f) can be de- 
composed in general as a "sum" of many symmetric ba- 
sis vectors |5(a;,a;)),the above observation means that 
the destructive phenomenon in two photon coincidence 
is caused by the symmetric components \S{u),u))) of the 
light field. Correspondingly the two wave packet is a 
"double sum" of 'i'^^d! over (uj.uj). As D = 0, every "i'^-a, 
approaches to zero and then \I' is also zero. Since the null 
phenomenon ^'|d=o = was contributed by every com- 
ponent '^ujCj\d=o = 0; we are led to the universal conclu- 
sion that there does not exist an effect of "interference 
between the two photons" . This conclusion holds for the 
separable case with a symmetric factorized distribution 
function f(io,Co) = namely, the two indepen- 

dent pulses have the same shape exactly. In this case 
this destructive interference phenomenon of two photon 
coincidence may occur due to the pairing mechanism just 
mentioned above. 



III. FREE INDUCTION DECAY MECHANISM 
FOR DESTRUCTIVE INTERFERENCE 

In the remaining part of this letter, we will use the 
free induction decay mechanism to explain the anti- 
exponential decaying behavior beside "dip" (or anti- 
peak) point. It is easy to see that each component 
\S(Q., 17)) in the symmetric two photon state contributes 
the coincidence counting rate Rc with the oscillating term 



cos 



Here, the term oscillating with T has been omitted in 
the long T limit. Apparently, in this sense , Rc is 
an oscillation function of D with a constant frequency 

Then, Rc is zero not only at D = , but 

also at D = Dfe = {k = ±1, ±2, ...). The oscillating 
term '^^(^j of different frequencies in the two photon wave 
packet can cancel each other beside the points D = 
and Dk, this will explain the "dip" experimental curve 
of non-oscillation in the two photon coincidence. Based 
on this conception of the free induction decay mechanism 
[15], the calculation of the coincidence counting rate Rc 
leads to the explicit result for the general symmetric in- 
put state 



r-oa jj 

RcOC dLi)du)\Ff {1 — cos[{uj — Q) —]} 

J-oo 



(9) 



Here, the integral domains of variables t{, uj and w 
have been extended to (— oo. +oo) for Gaussian type dis- 
tributions f{LO,ui). It is easy to see that each oscillating 
term of variable D in Rc has a frequency {ui — Hi). If 
\F\'^ is not very steep at a certain point (wq, wo),the dou- 
ble sum of cos [(w — w) ^] with weight \F{co,lj)\'^ will 
become a decaying function of variable D with the maxi- 
mum at _D = and then Rc will become an anti-decaying 
function. This is a typical optical free induction decay 
phenomenon for the second order coherence function. 

For an illustration of the two photon free induction 
decay, let's calculate a separable case with F{uj,lo) = 
F{uj)F{uj) where F{u}) is of Gaussian type, i.e., F{co) = 

^"2 "2°^ ]• It is straightforward to calculate Rc , 



exp 

obtaining the result Rc oc 



From this 



calculation we observe that, even for a factorizable two- 
photon input state, there still occurs the destructive phe- 
nomenon Rc\d=o = with the anti-exponential Gaus- 
sian decay as \D\ increases if the two photons are cre- 
ated at the same time and their shapes are strictly the 
same i.e ., f{ui) = f'{u)). It is remarkable that this phe- 
nomenon i.e. the interference of two independent photon 
has been analized theoretically and realized experimen- 
tally recently [23,24]. 

To investigate the possibility of testing the present the- 
oretical predictions, we have consider the practical cases 
where the condition G{uj,uj) = is not strictly satisfied. 
Obviously, just at = 0, the counting rate of two photon 
coincidence 



J — C 



dLodCb\G{uj^Cb)\^ 



(10) 



(8) 



measures the average extent of asymmetry of the input 
two photon states, which depends on the difference be- 
tween F{u},Lu) and F(Gj,uj). On the other hand, around 
the point £> = 0, 
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Rc = J |*lpdT{dT^+ j |^2pdT{dT^ (11) 

+2Re J ■^l^2dT[dT^ 

It is easy to prove that, as increases, the first term in 
Rc will increase while the second term keeps a constant 
independent of D. However, the presence of the third 
term may suppress the increasing tendency of i?c as \D\ 
increases. In the case that |G(w, (I>)|is very small and 
1^2! << 1^1 17 the shape of Rc in the neighbor around 
Z? = is just similar to that for the ideal case with 
G{u,u) = 0. 

In Fig3, we give the curve of "dip" for the rate of 
two photon coincidence in a typical non-symmetric case 
that G{lu,u}) ^ and F{uj,u}) — Fi (0^)^2(0)) where 

i^iH = exp(-i^) and F2{Cj) = exp(-(^). The 
asymmetry parameter y is defined by y = 1 — j|. It de- 
scribes the relative difference between ft and il. In this 

2 r-,2 2 2 

case we have i?c oc 1 — exp(— ^2^2 ^^^)- ^^^V 

to see that, as y increases, the value of Rc\d=q increases 
and " the Rc— curve" becomes more and more flat. 
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FIG. 4. The curve of Rc = r as a function of D when 

F{lo,lo) — e 2<T^ e e ^'^ . Here we define parame- 

ter y = Ta. The unit of 73 is f . 

So far we have only discussed the ideal situation with 
50 — 50 beamsplitter where the reflectivity (i?) and trans- 
missivity (T) of the beam splitter are all 1/2. In the 
general case Ry^T, we have 

E+ {h ) = iVREt (ri ) + VtE+ {r[ ) (12) 

and the corresponding two-photon 

wave packet^(ti, t2).In this general case, the necessary- 
sufficient condition for Rc\d=o — can never be satis- 
fied unless i? = T = 1/2. In the experiment [9] where 
G {uj,uj) = and i? 7^ T, the effective wave packet \E' can 
also be decomposed into two parts:\E'i cx R dujdu}{...) 
x(l-e"*(""'^)^),*2 oc (T-i?). Apparently, when £) = 0, 
^1=0 and Rc\d=q — |^2p- Therefore, ^'2 characterizes 
the basic difference between R and T. If |y — 1| is very 
small, ^'2 can be considered as a perturbation for '^i. 



FIG. 3. The curve of i?c = r as a function of D when jy TWO PHOTON INTERFERENCE OF TYPE 

F{lo,lo) = e 2^^ e 2^^ . Here we assume — = 10. The 

unit of D is " II 



Another more interesting illustration is that F(lu, Co) — 
Fi{uj)Fi{u})e^'^'^'^ . This distribution function means that 
the two-photon wave packet can be factorized as the 
product of two signal photon wave packets of the same 
shape, one (in mode A)oi which was created earlier than 
another (in mode B) with time T. In this case, it is easy 
to prove that Rc{T, D) = Rc{0, D + cT) . In Fig.3 . 



The above analysis can be generalized to the case of 
photon's polarization. For instance, we can consider the 
two-photon interference demonstrated in Fig. 5. 




FIG. 5. Schematic of the two-photon interferometer with 
polarization considered 
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In this theoretical scheme of two-photon interference 
experiment, a photon pair state with entanglement with 
polarization degrees of freedom can be written as 



where 



I*) = / 

Jo 



(kodivf {lo, lv) X 



(13) 



The photons arrive detector Di and D2 after mixed by 
the beam splitter BS. Here [uo) is the creation op- 
erator of photon of frequency w, the wave vector is along 
direction Oa (in optical path mode A) with horizontal po- 
larization (along direction ch)- The notations a^y (S), 
'^iy '^^'^ '^■BH ^'''^ defined similarly as 
Two polarization analyzers Ai and A2 along the direc- 
tions Ai and A2 are put in front of the detectors. 

It is easy to see that the rate R,. of coincident detec- 
tion is proportional to the time integral of the square of 
effective wave function ■^'s mode, i.e. 



where 



Rcoc J dtidt2\tp{ti,t2)\^ 



^P{h,t2) = {0\E+{h)E+{t2)\^). 



(14) 



Here E+ (U) = A,-E+ (i,),and E+ {U) = ch E+ (t,) + ey 
Ey {ti){ i = 1,2) are the positive frequency parts of light 
field operator at the position of the detector. Considering 
the transforming character of beam splitter, we have 



Et{U) = ^[iEt{T,) + E+{T[) 
Eiit2) = ^[^{r!,)+iEl{r2) 



(15) 



Here n, r{, T2, T2 are just defined as in section 2 and 

Ea{b) = en + E+^^y are the positive fre- 

quence part of the light field operator at point a and b 
^^'^ Kib)H{v) (*) = / d<^(^A(B)H(v) exp [-iujt]. The spec- 
tral density function g (u) is neglected here without loss 
of gcncrity. 

Submitting 15 to 14, we can get the explicit expression 
for effective wave function ip of the light field's : 



1 



dwdujf X 



A-i ■ ev X A2 ■ bh - Ax ■ en ^ A2 ■ ey 



X (16) 



/ 



where D = li — si = c{ti — t[) = c (tj — T2) is the op- 
tical path difference disscussed before. Then can be 
written as 



(17) 



Ai-ev X A2-eH - Ai-bh X A2- ey 

(18) 

is dependent of the direction of polarization analyser and 
r]{D) = j dT[dT2\ j dwdwf x 

g-i(air{-|-wr^) _|_ g-j[a;r{-|-wr^-|-(w-w)-|] |2 ^^9) 

varies as the optical path difference D is changed. 

It is easy to prove that, when the frequency distribu- 
tion function / is symetric, i.e. / (w, w) = / (w, w) , we 
have ^ |D=oand 



dD 



dD^ 



|d=ooc — j dijjduj l/l^ (w — w)^ < 0. 



(20) 



So we know that in this case Rc has its maximal value 
when D = Q. This is differenct to the case in which 
the photon polarization is not considered and i?c has its 
minimal value when D = 0. In fact, this difference is due 
to the —1 factor in light field's initial state 13 i.e. if the 
initial state is 

dwdwf X --j= [a^H (w) a^y (w) + a\y (w) a^^ i^)] |0) ' 

(21) 

Rc will also has its minimal value when D = 0. 

It is notable that when the photon's polarization is 
considered, / {w,^) = f (w, w) is only the sufficient con- 
dition of that Rc has its maximal value (or minimal value) 
when D = but not the necessary condition. In fact, it 
is difficult to obtain its necessary condition. 



V. CONCLUDING REMARKS 

To summarize, we have theoretically investigated the 
two photon interference phenomenon for a general situ- 
ation with an arbitrary two photon state and given the 
necessary-sufficient condition for the destructive interfer- 
ence of coincidence counting in terms of the symmetric 
pairing of photons in the light pulses, i. e. , G(a;, w) = 0. 
This implies that, even in the separable case f{uj,u}) = 
fi{'^)h{'^), if /i = /2,a curve of "dip" can still be ob- 
served in the rate of two photon coincidence. Our investi- 
gation predicts the destructive interference phenomenon 
to occur even for certain cases with separable input two 
photon state. However, this does not imply "the interfer- 
ence between two independent photons" since the essence 
leading to destructive interference is that the required 
two photon state consists of the inseparable symmetric 
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componcnt|S'(ti', cj). Since for the the symmetric distribu- 
tion function all components of the light field state have 
the entangled forms |5(ci;,ci;),we understand the "dip" in 
the destructive coincidence curve according to the free in- 
duction decay mechanism from the dispersion of the two- 
photon frequency {uj,lj). In fact, though each component 
contributes Rc with an oscillating term with respect to 
D, in the integral of variables w and u with Gaussian type 
distribution / , these oscillating terms cancel one another 
and thus lead to a (anti-) Gaussian decaying factor. We 
also considered the effect of photon's polarization. We 
found that, when the effective wave function is symetric 
i.e. f {uj,Lu) = f (Cj,lo), the coincidence counting rate 
may have its maximal or minimal value when the path 
difference is zero. 
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